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1. Introduction 

Since long it is suspected that SU(A^) gauge theories can be regarded as the low energy ef¬ 
fective theory of N D-branes in specific parameter regimes. In this way, dimensionally reduced 
large-A^ supersymmetric Yang-Mills (SYM) gauge theories might provide a nonperturbative for¬ 
mulation of the string/M-theory which describes the dynamics of the D-branes. In particular, the 
connection between black p-branes and SYM gauge theories in <i = (p -|- 1) dimensions allows to 
study black hole thermodynamics through the corresponding strongly coupled gauge theory. Apart 
from providing tests of the gauge/gravity duality, SYM quantum mechanics is interesting per se. 
There are various intriguing expectations concerning the behaviour of the theory in specific secfors 
of fixed fermion number [1,2]. In particular, while certain canonical secfors only have a discrete 
energy specfrum, some secfors allow for fhe discrefe specfrum fo be immersed in a confinuous one 
reaching down fo zero energy. These confinuous specfra can presumably be associafed with the 
so-called flat directions of the potential which are present in the classical theory, but may or may 
not survive the quantisation of the theory. 

In these proceedings we report on our ongoing effort to perform nonperturbative calculations 
in SYM quantum mechanics with gauge group SU(A). Here we concentrate on = 4 SYM 
quantum mechanics which is obtained from dimensional reduction of the jV = 1 SYM gauge 
theory in d = A dimensions by compactifying the thr'ee spatial dimensions. To define fhe fheory 
we employ fhe lattice regularisafion in Euclidean lime proposed in [3, 4]. The bosonisalion of fhe 
fheory on fhe lattice by means of fhe fermion loop formulation [5, 6, 7] decomposes the fermion 
contributions into fermion sectors with fixed fermion number. Our recenf progress in underslanding 
fhe algebraic sfrucfure of the fermion loop formulation allows the explicit construction of transfer 
matrices [8]. The transfer matrices in turn provide the starting point for the construction of local 
fermion update algorithms which allow to directly simulate fixed canonical sectors of fhe fheory. 

In fhe following we summarise fhe derivation of the transfer matrices for generic fermion num¬ 
ber sectors, recapitulating our results from [8], and show first results for some simple observables 
such as the Polyakov loop, the moduli of the bosonic fields and the fermion action from simulations 
in fixed canonical secfors of fhe theory. 

2. Lattice regularisation and canonical sectors 

We directly start with jV = 4 SYM quantum mechanics obtained from dimensional reduction 
of JC = 1 SYM ind = A dimensions down to r/ = 1 dimension. The action in Euclidean time can 
be written as 



(2.1) 


where Df = dt — i[A{t),-] denotes the covariant derivative using the time component of the SU(A) 
gauge field A(t), while fhe spafial componenfs become fhe bosonic fields Y,(t) wifh /= 1,2,3. The 
anficommufing 2-componenl complex fermion fields t/r(t) inferacf wifh fhe bosonic fields Y,- 
fhrough a Yukawa-fype inferacfion involving fhe fhree Pauli-mafrices a,. We nofe fhaf all fhe fields 
are in fhe adjoin! represenfafion of SU(A). The discrefisafion of fhe Eagrangian on a fime laffice 
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with Lf points is straightforward and yields for the bosonic part 

= 4 L Tr j DMt)DtXiit) - ^ [Xiit),Xj{t)A (2.2) 

^ ?=o t ^ J 

where DtXi{t) = U{t)Xi{t + l)f/^(t) — X,(t) is the covariant forward derivative and U{t) is an ele¬ 
ment of the gauge group SU(A^). For the fermionic part of the action we introduce a Wilson term 
with Wilson parameter r = 1 in order to avoid fermion doubling. With this choice of r the massless 
Wilson Dirac operator in one dimension involves just the forward derivative and one obtains 

I L,-\ 

= [Xi{t)Mt)]} (2.3) 

S t=o 

for the fermion action. More specifically, we have 


Sf 


1 

2g2 


I 


t=0 


-Ta mfp (^ + 1) + Wa it)Vp (0 


YS!p,a[U,Xi-,n]Y (2.4) 


where W{t) denote the real adjoint gauge link matrices 


W5(0 =2(ao)«,3®Tr{rt/(0r''^(0^}. 


(2.5) 


Note that we have introduced a chemical potential /r in the standard way [9]. The subscripts p,a 
for the Dirac matrix ^ denote periodic or antiperiodic boundary conditions, respectively, for the 
fermionic fields. The Yukawa interaction matrices <I>(t) are x with = 2{N^ — 1) and 
read 

Kpit) = (ao)a/3®5“^-2(a0a/3®Tr{r[X,.(O,n}. (2.6) 

Two remarks are in order. Firstly, all supersymmetry breaking terms apart from the ones 
introduced by the discretisation in eq. (2.2) and (2.4) are forbidden by the gauge symmetry. Hence, 
supersymmetry is expected to be automatically restored in the continuum limit without any fine 
tuning [3, 4]. Secondly, the Wilson term breaks the time reversal symmetry, or equivalently the 
particle-hole exchange symmetry. This reflects itself in the fact that the action in eq. (2.4) only 
allows forward propagating fermions. As a consequence, the exchange symmetry between the 
related fermion sectors with and — tif becomes exact only in the continuum limit. 

Let us now derive exact expressions for the fermionic contributions to the partition function 
of the theory for a given fixed gauge and boson field background - the canonical determinants. In 
quantum mechanics the lattice regulated determinant of the Dirac matrix can readily be calculated, 
and one obtains 

L,-l 

det^p,a[f/,Z;;At] =det[5^=Fe+^^'] with ST = {t). (2.7) 

r=0 

This essentially corresponds to the dimensionally reduced determinant for Wilson fermions derived 
in [10, 11] except that here the dimensional reduction is from the full matrix to a x 
’flavour’ matrix. It is now easy to get the canonical determinants from the fugacity expansion 

2 (N^-\) 

dtt%p[U,Xi-ii\= ^ (2.8) 

«/■—0 
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which identifies the canonical determinants as the coefficients of the characteristic polynomial. 
These coefficients can be expressed in terms of the elementary symmetric functions Sk of order k 
of the eigenvalues {t,-, / = 1, • • •, of 

k 

Sk{^) = Sk{{i:i})= X ( 2 . 9 ) 

and one eventually obtains 

det%[f/,X,] = . (2.10) 

So the crucial object for the calculation of the canonical determinants is the product ^ of the 
matrices <I>(f) and W{t) which in fact is a product of transfer matrices [8], as was suspected already 
in [10]. 


3. Transfer matrices for the canonical sectors 


The explicit construction of the fermion transfer matrices for each fermion sector is most easily 
done via the fermion loop formulation [8] which in essence is an exact (hopping) expansion of the 
fermionic Boltzmann factor to all orders. In this formulation, the contributions of the fermions 
to the partition function are obtained by summing over all possible closed oriented fermion loops 
which are forward propagating in time for any given gauge and boson field background. The 
loop configuration space naturally separates into subspaces characterised by the number of forward 
propagating fermions nf. The transfer matrix elements are explicitly given in terms of the cofactors 
C(d>) and the complementary minors M{W), 

= C^^(<I>) = (-l)^('^>«)det<I>^H, (3.1) 

( 27)^6 ^ 

where A,B are sets of indices A,B C {1,... of order Uf and p{A,B) — HieA * + L;eB j- 

denotes the matrix obtained from <I> by deleting the rows with indices from B and the columns 
with indices from A, while denotes the matrix obtained from W by picking only the rows with 

indices from A and columns with indices from B. The size of the transfer matrices is given by the 
number of such sets for a given n j and corresponds to the number of forward propagating fermion 
states A^states = !/(uj— «/•)! •«/!. The fermion contribution to the partition function in sector 

nf is then simply given by 


det%[[/,X;] =Tr 


1177 ( 0 - 77(0 


(=0 


(3.3) 


and one can use the Cauchy-Binet formula and some further algebra [8] to show that 


nr7(0-r7(0 

t=Q 


(_l)riA5)det,^X^ = C^^(,^), 


(3.4) 


4 







Canonical simulations ofSYM quantum mechanics 


Urs Wenger 


CL 




0.8 


0.7 


o Dpi 6 

o n,=15 

A n,= 1 
□ n,= 0 



0.6 


0.5 


0.4 



0.3L 


_ 1 _ 

0.5 


1 




P 




f 11 


o 

n,=16 


n,=15 

< 

n,=14 

> 

n,=13 

A 

n,=12 

O 

n,=11 




i|ii« 






Figure 1: The modulus P of the Polyakov loop as a function of j3 for SU(3) on a L, = 5 lattice in various 
canonical sectors. 


hence the canonical determinant is simply given by the sum over the principal minors of order n f 
of =^7 denoted by En^{£Y), 

det [U,Xi] = £ det = Enf{^). (3.5) 

B 

Finally, it is easy to show that E„f{^) = which proves the equivalence between the 

representation using the transfer matrices and the one in eq. (2.10). 

Some remarks are in order. Firstly, we note that the matrix ^ describes the dimensionally 
reduced effective action for the Polyakov loop coupled to the bosonic fields X,. Secondly, our 
result for the canonical determinants in principle allows for local fermion update algorithms, but in 
practice only the sectors with u/ = 0 and «/ = can be implemented straightforwardly, while 
in other sectors algorithms along the lines in [12] can be employed. Thirdly, the construction of the 
transfer matrices and the calculation of the canonical determinants in terms of those is applicable 
to QCD, since the algebraic structures of the theories are the same. 


4. Canonical simulations 


Here we present our first results from simulations of the system employing the gauge group 
SU(X) with N = 3 and = 16 directly in the various canonical sectors ^ First we note that 
the canonical determinants are real because the eigenvalues T; of ^ are real or come in complex 
conjugate pairs. Furthermore, for the sectors nf = 0 and «/ = (quenched) one can prove that 
the canonical determinants are positive. In these two sectors we update the bosonic degrees of 
freedom using a local Metropolis algorithm. In the other sectors we use Metropolis updates based 
on eq. (3.5) and currently simulate only in the configuration space with positive determinants. 

In the following we show results on a lattice with temporal extent L, = 5. We measure the 
moduli of the Polyakov loop and the scalar field defined by 


P = 


TrnC(0 

t 


= \X\'^ =X^X^. 


'Note also the recent effort using the RHMC algorithm for SU(2) in [13]. 


(4.1) 
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Figure 2: The square of the modulus of the scalar field as a function of j3 for SU(3) on a L, = 5 lattice in 
various canonical sectors. 


We note that in some sectors the simulations become unstable and grows without bound. We 
believe that this is because the flat directions become unstable due to lattice artefacts, and we expect 
the behaviour to disappear towards the continuum limit. In the left panel of Figure 1 we show P as 
a function of the temporal extent of the system parametrised by j3 for the sectors n/ = 0 and 16, and 
for rif = I and 15. Each pair should be degenerate in the continuum and we see that this is indeed 
the case for jS < 0.6, while for larger values of j3 lattice artefacts lift the degeneracy. The (physical) 
differences between the sectors is illustrated in the right panel of Figure 1 where we show P in the 
sectors tif = 16 down to 11. 

Next, in Figure 2 we show the square of the modulus of the scalar field as a function of jS 
for the same combinations of sectors as in Figure 1 . Again we find that the degeneracy between the 
Jif = 0 and 16, «/ = 1 and 15, and so on, is lifted by lattice artefacts towards large values of j3. 

Finally, in Figure 3 we show the fermionic action Sf = (In det &nf)nf as a function of jS in 
various canonical sectors. We find that the degeneracy of this observable between the mirror sectors 
becomes better and better towards jS —>• 0, suggesting that reweighting between the min'or sectors 
could become feasible in that limit, or more generally towards the continuum limit. 
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Figure 3: The fermionic action as a function of j3 for SU(3) on a Lf = 5 lattice in various canonical sectors. 
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5. Summary and outlook 

In this contribution we summarise the derivation of explicit transfer matrices for ^ = 4 SYM 
quantum mechanics with generic gauge group SU(A^) discretised on a time lattice. The transfer 
matrices are defined separately in each canonical sector with fixed fermion number n / and form 
fhe basis for canonical simulafions of fhe fheory. One caveaf is fhaf in fhose seclors where fhe 
canonical determinants are not positive definite, the local Metropolis algorithm is currently not 
very efficient and only samples configurations with positive determinants. 

Several paths are now open for further investigation. From an algorithmic viewpoint, it is in¬ 
teresting to examine the systematics of reweighting ensembles of configurations from one fermion 
sector to another, or from simulations at finite (imaginary) chemical potential. Concerning the 
physics of the model, it is interesting to calculate correlation functions and energy spectra in the 
various canonical sectors. The investigation of the phase transition in the large-M limit of the 
= 16 SYM quantum mechanics is most useful for a further understanding of the thermodynam¬ 
ics of certain black holes. The results in these proceedings are a first step towards these goals. 
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